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Abstract A review is given for cyclotron resonant interactions in space plasmas. After 
giving a simple formulation for the test particle approach, illustrative examples for res- 
onant interactions are given. It is shown that for obliquely propagating whistler waves, 
not only fundamental cyclotron resonance, but also other resonances, such as transit- 
time resonance, anomalous cyclotron resonance, higher-harmonic cyclotron resonance, 
and even subharmonic resonance can come into play. A few recent topics of cyclotron 
resonant interactions, such as electron injection in shocks, cyclotron resonant heating 
of solar wind heavy ions, and relativistic modifications, are also reviewed. 
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1 Introduction 

It is generally believed that for efficient nonthermal particle acceleration in various 
explosive cosmic events strong plasma turbulence play the essential role. For non- 
thermal particles to emerge from the pools of majority thermal particle populations, 
some selection mechanism(s) should work in turbulent environments. The most likely 
selection mechanism is due to resonance interaction which assures the efficient en- 
ergy/momentum transfer from the turbulence to nonthermal particles. In this review, 
we will discuss the cyclotron (and transit-time) resonance interactions which would be 
viable in various space plasma environments. 
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2 Test particle approach 

In usual space-plasma environments, the electromagnetic fields have the background 
components (Eo, Bo) and overlapping multi-wave components, (EJ^, BJ^) where # = 
1, 2, ... are the mode ID numbers. The generalized Ohm's law gives the relation between 
Eo and Bo as 

E = -(l/c)V wfc x B , 

where we write the bulk velocity of the plasma as ~Vbulk an d neglect higher order terms 
caused by the spatial variations of Bo and electron thermal pressure. If we make the 
coordinate transformation to the rest frame of the bulk plasma, Eo vanishes. For a 
nonrelativistic electron, the equation of motion i£] 



where Q C e is the electron cyclotron frequency eBo/m e c. JT]) can be regarded as an oscil- 
lator equation of the eigen-frequency Q C e with the forcing term on the right hand side. 
When we integrate ([l]), we evaluate EJ^ and B^ at the electron position r e (t) which 
is affected by the existence of the waves. Therefore the system of equations JT]) and 
([2]) have nonlinearity, from which mode coupling effects among the wave components 
can appear. However, if the amplitudes of waves are small enough, we can separately 
integrate (JTJ) for each wave, 

q± + Qce W e x|l = -^{Ef + i Ve x B#| (3) 
at |B | m e I c J 

and then superpose the results afterward (the standard procedure in the quasi-linear 
treatment). From here we choose a single wave mode and drop the superfix #0- 

For simplicity we consider electromagnetic waves in a cold uniform plasma con- 
sisting of electrons and ions (protons) of equal density no- We define the frequency 
and wave number of the waves as to and k, respectively. If not otherwise noted, the x 
axis is taken along the background magnetic field Bo, and the z axis is so chosen that 
the wave vector k is in the x — z plane, k = (k x , 0, k z ) = k(cos0, 0,sin#). Linear 
dispersion property is summarized in Appendix A. We neglect the back reaction of 
electrons toward the electromagnetic field (test-particle approach). 



1 Here we have chosen the electron motion. Of course, we can equally take the ion motion 
with a suitable modification of the choice of wave modes in the following consideration of 
resonant interactions. The relativistic effect will be considered in £15.41 and £15.51 

2 Note, however, that there still remains self nonlinearity in J3j for the mode #, where 
E* and B* are evaluated at the particle position r e (t), which is affected by and B* 
themselves. This nonlinearity is the origin of subharmonic resonances to be discussed in i|4.5l 
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Fig. 1 The result of test-particle calculation for a parallel propagating wave with oj = 5J7 c i = 
2.72 X 10~ 3 J? ce , k = k x = 2.0iAQ ci /V A , and \B W \/B = 10~ 4 . The variance of /j,, the cosine 
of the pitch angle, is plotted against the average velocity (V ex ), where the horizontal axis is 
normalized by the Alfven velocity Va- The peak at (Vex) /Va at ~ —892.2 shows the effect of 
the fundamental cyclotron resonance condition. 



3 Parallel propagating waves 

The first case is for a parallel propagating whistler wave with oj = 5f2 c i = 2.72 x 
10 i? C e, k = k x = 2.044J? c i/V^4, and the relative wave amplitude |Bu,|/i?o = 10~ 4 . 
We inject 3.5 x 10 4 electrons with the initial conditions of 

r e = (0, 0, 0) and V e = (V e x.oi V e ±cosa, V e ± since), (4) 

where (V ex0 , a) are distributed uniformly in the two dimensional space, (—GOOOVa, 
6000Va) x (0, 2tt), and V e< ± is fixed at IOOOVa- We integrate @ and © for a long 
enough time T, and then evaluate the averages of Vex, M (cosine of the pitch angle), 
and 8/j,8n (the variance of /x), 

(Vex) = i / Vex(t)dt 

1 Jo 

In Figure Q] (SfiSfj,) is plotted against (Vex), where we see a peak of (SfiS^j.) at (Vex) = 
— 892.21^ = Vr, where Vr is the resonance velocity satisfying 



UJ — k x Vji = nflce 



(6) 
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Fig. 2 The same as Figure [T] except that it is for a oblique propagating wave with uj = 
5f2 ci = 2.72 x l(r 3 i? ce , k = 2Ann ci /V A , k x = 1.709J2 ci /VA, 8 = 45°, and \B W \/B = 
10~ 4 . Integers given above the peaks (—5, —4, , 5) are harmonic numbers in the resonance 
condition, u — k x Vn = nQ c£ . Half integers given below the curve (—5/2, —3/2, ...,5/2) are for 
the subharmonic resonance (see text). 



with n = +1 (the fundamental cyclotron resonance condition). This is the most natu- 
ral resonance, in which right-hand rotating electrons interact with right-hand polarized 
whistler waves. Historically-famous Kennel-Petscheck theory 1 Kennel and Petschekl (|l966l )) 
for the loss-cone instability of whistler waves by t he trapped electrons in the radiation 
belt mainly focused on this resonance. Recently. lAmano and Hoshinol ( 20ld ) applied 
this resonance condition to the electron injection process at shocks (^SHJ. 

Note that electrons in a monochromatic wave are 'trapped' around Vr and their 
pitch angle diffusion does not occur. In Figure 1 there is a dip at {Vex) = Vrj which 
is the manifestation of the trapping effect (Appendix B). If, on the other hand, we 
superpose multi waves, the interaction between the electron and one particular wave 
component is limited within a finite coherent time r, giving a contribution (5(i5[J,}/t to 
the pitch-angle diffusion coefficient. The task of quasilinear theory from this point 
of view is to evaluate the coherent time r for a given ensemble of wave compo- 
nents. Since the literature on quas ilinear theory is too vas t , we s h ould limit the ref- 
erence to only a few articles (e.g., 



Kennel a nd Engelmann Jl966h: Hall and Sturrockl 



19671) :lHasselma,nn and Wibberenz ( 19681) :lKulsrud and Ferraril dl97lh:IJokir>iJ dl97lR 



Lee and Lerchd h97^;ISkilhnd (Il97^ and Isenberd 

(20021); ISchlickeiserl (|2002h : IPetrosian and Liul (|2004l ); IShalchil l|2009l )). 



4 Oblique propagating waves 

Figure[2]is for an oblique propagating whistler wave with uj = 5J? c i = 2.72 x 10~ 3 i? C e, 
k = 2. 417 fi c i/VA, k x = 1.709O c i/VA, 9 ~ i"/4, and the same relative wave amplitude 
|Buj \/Bq — 10 -4 as in We inject 5. Ox 10 4 electrons with the same initial condition as 
0. Peaks labeled with integer numbers (n = —5, —4, ,5) correspond to the harmonic 
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Fig. 3 Schematic illustration of the three resonance conditions, U3 — k x Vn = nQ ce , for 
(n = 0, ±1) for an electron with a whistler wave (red circle). 



number n in the resonance condition (|6} . It is noted that there sometimes appear dips 
at the resonant velocities. Again these dips are manifestation of the trapping effect 
already seen in Figure [1] The resonance conditions for n = 0, ±1 are illustrated in 
Figure [3] 

4.1 Fundamental cyclotron resonance with n — +1 

The highest peak with n = +1 in Figure [5] is for the fundamental cyclotron resonance 
condition, 

w - kxV R = +f2ce (7) 
which is already discussed in ij3]for the purely parallel propagating case. 

4.2 Transit-time resonance with n — 

The second highest peak in Figure [5] is for the Landau resonance or transit-time reso- 
nance condition, 

w - k x V R = (8) 

To see the origin of this resonance, we return to the x component of the equation of 
motion (O, 

— 77— = { E w .x H (Ve,yB w . z — V e ,zB w .y)\ (9) 

where the second term in the parenthesis on the right hand side of ([2} is the (V e x H w )x 
term. From the numerical inspection of ([9|, we find that the relative contribution of 
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the E W: x term is less than 1% of the (V e x B K )i term. (The dependence of {5)x5fj) 
on (Vex) is essentially unchanged, if E w , x is artificially set zero.) The (V x B) x term 
describes the mirror force exerted by the variation of t he magnetic field magni t ude, a nd 
this resonance is called 'transit-time' resonance (e.g.. lSchlickeiser and Miller! i|l998l )). 



4.3 Anomalous cyclotron resonance with n = — 1 

The resonance corresponding to the peak with n — — 1 in Figure [2] 

w - k x V R = -n ce (10) 

is hi storically c alled ' anoma l ous cy clotron resonance' (or 'anomalous Doppler effect': 

to 



is historically called anomalous cyclotron resonance ( 
e.g. JCinzburgl (|l960h ; iBricel (fl964)1. (flUTl is equivalent 



f - V * = ~T L (<0) 

fax fax 

which means that the wave is overtaken by the electron as Vr exceeds the wave phase 
speed ui/kx (here we set k x > 0). Therefore, there is a reversal of the wave polarization 
from the background plasma frame to the electron comoving frame. The term 'anoma- 
lous' stems from this reversal of polarization. It is noted that an oblique propagating 
whistler wave is elliptically polarized and consists of right-hand polarized (R) and left- 
hand polarized (L) subcomponents. When the resonance condition (|10|) applies, there 
are dual polarization reversals, namely R and L subcomponents in the background 
plasma frame respectively become left-hand and right-hand polarized in the electron 
comoving frame. What the electron efficiently exchanges energy and moment um at the 
n = — 1 condition is the L subcomponent in the background plasma framj^l iLevinsonl 
( 19921 ) utilized this resonance condition to excite oblique whistler waves by electron 



beams injected from the shock front. 



4.4 Cyclotron higher-harmonic resonance at n with |n| > 2 



Cyclotron higher-harmonic resonance phenomena for o blique and perpendicular prop- 
agating waves are well known (see, e.g., Chapter 10 of IStbJ (1992)). Its application to 
Bernstein wave observed in the Jovian magnetosphere is one of the classical examples 
( Barbosa and Kurtbl (|l980h i. 



4.5 Cyclotron sub-harmonic resonance at half integers 

In Figure [2] there are some minor structures at subharmonic numbers, 

n = -5/2, -3/2, ..,3/2, 5/2, 

where the dip structures (trapping effects) are seen to dominate over the small peak 
structures. These are manifestations of cyclotron sub-harmonic resonance interaction, 

3 From this consideration it is naturally understood why there is no anomalous resonance 
effect for a parallel propagating whistler wave (Figure 1): An electron satisfying the n = — 1 
resonance condition feels the whistler wave left-hand polarized in its own comoving frame and 
does not exchange energy and momentum with the wave efficiently. 
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Fig. 4 Angular dependence of (SfiSfi) for whistler waves at three different resonance condi- 
tions (n = 0, ±1) is shown for 9 = 0° ~ 85°. The wave frequencies and relative amplitudes are 
fixed at u = 50 ci = 2.72 X 10 _3 i? ce , and \~B w \/Bq = 10 -4 , respectively. 



wh ich was first discussed b y Smirnov and Frank-Kamenetskiil ( 1968). and then followed 
by iTerasawa and Nambul l|l989h for perpendicular propagating magnetosonic waves. 
There are r e newed interests on this partic u lar re s onance for obli q ue Al fven waves (e.g., 
IChen et al.1 (|200ll ); IWhite. Chen and Linl l|2002l ); iLu and Cheiil (|2009h V The origin of 
this resonance is the self nonlinearity noted in the footnote^ of ij2] 



4.6 Angular dependence of {S/iSfi} 



As shown in Figure 2, there are many possible resonance modes for oblique propagating 
waves. It is of interest to see how the relative importance of these mode depends on 
the propagation angle 9. Figure 4 shows the dependence of (<5p<5/i) on the propagation 
angle 6 of whistler waves at three different resonance conditions (n — 0, ±1). In the 
quasi-parallel propagation regime (9 <~ 45°), the fundamental resonance of n — +1 
dominates over others. On the other hand, in the quasi-perpendicular regime (9 >~ 
45°), the transit-time resonance (n = 0) and the anomalous cyclotron resonance (n = 
— 1) become progressively more important as 9 — > 90°, and exceed the fundamental 
resonance above 9 ~ 60° — 70°. 
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Fig. 5 From lAmano and Hoahind ll201df) . Schematic dispersion diagram for whistler and ion- 
cyclotron waves in an electron-ion plasma. Positive (negative) frequency corresponds to the 
right-hand (left-hand) polarization. The cyclotron resonance condition J7J is also shown. Waves 
in the shaded regions are strongly damped by the cyclotron damping of thermal plasma par- 
ticles (electrons and ions). 



5 Recent topics of resonant particle interactions 

5.1 A critical Mach number for electron injection in shocks 

It has been known that the e lectron component of galactic cosmic rays is about 1% of 
the proton component (e.g.. ISchlickeiserT l(2002)), While the origin of this percentage 
is not well understood, it is generally interpreted as the indication of smaller injec- 
tion/acceleration efficiency for electrons than ions. For the case of ions, a part of their 
suprathermal component is reflected by shock fronts to form a beam. The beam ions 
then excite MHD waves through the ion-bea m-cyclotron resonance con dition and are 
subject to the pitch angle scattering process ( Winske and Lerovl fl984; ) ) , which is the 



first step to the diffusive shock acceleration process for the ions. If we apply the same 
scenario to electrons ; difficulty arises at the level of the wave excitation. Figure [5] illus- 
trate this difficulty l Amano and Hoshinol l^Old )): An oblique solid line running from 



the upper-left is the fundamental cyclotron resonance condition J7J for the electrons 
with the beam velocity, V e beam: 



(12) 
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Fig. 6 From lOka" et alJ l l2006f ). (a) Observed spectral indices F of suprathermal electrons 
at the bow shock transition region are plotted onto the two dimensional map of observed 
parameters (Ma, 8g n ). F is color-coded from orange (hard) to dark blue (soft). A solid black 
curve shows the position of the critical Mach number Ma, crit where is set 1. (b) Observed 
spectral indices F are plotted against the Mach number normalized by Ma crit where observed 
0Bn is used to calculate Ma, crit- 



The interaction point 'A' between the beam electrons and the whistler wave is not 
available for the wave excitation since the beam electrons with purely parallel velocity 
absorb whistler wave energy at this point. The interaction point 'B' is for the anomalous 
cyclotron resonance condition between the beam electrons and ion cyclotron wave 
(left-hand polarized in the plasma rest frame). However, in the interplanetary and 
interstellar plasmas of typical temperature, there are a plenty of thermal ions which 
damp the ion cyclotron wave via fundamental cyclotron resonance interaction (shaded 
region ar ound 'B') so that effic i ent pi tch angle scattering of electrons is prohibited. 
What lAmano and Hoshinol ( 20ld ) proposed is to invoke the loss cone distribution 



of beam electrons, which is the natural consequence of the electron reflection process 
at the shock front (shock drift acceleration process) and can provide the free energy 
to the whistler wave excitation at the interaction point 'A' if the damping effect by 
thermal electrons (shaded region around the u) = O ce line, namely —k x v e u l + fi C e < 
ui < fe x i! e tft + Sice with the thermal velocity v e tfr of background electrons) is avoided. 
These considerations lead to the criterion for the shock Alfven Mach number M A for 
the efficient electron acceleration, 



M A >M A>critS ^^J^0 e (13) 
2 y m e 

which agrees with the observational criterion for electron acceleration obtained at the 
earth's bow shock, where the hard spectral indices of suprathermal electrons (F < 3.5) 
are observed p referentially for h igh Mach number bow shock satisfying M A > M A cr j t 
(Figure |5] from lOka et all (|2006h ) . 



10 




10 20 10 40 

t(O) 

Fig. 7 Figures from Kuramitsu and Krasnoselskikh (2005a). Trajectories of ions satisfying a 
forced resonance condition with an electromagnetic wave under an external electrostatic field. 
Resonant (solid lines) and near-resonant (dashed lines) ions in (a) v x — x and (b) v z — v y 
phase spaces. A dotted line in (a) denotes the spatial profile of the electrostatic field, (c) Time 
evolution of v v of the resonant ions (solid lines) and E y (dotted line) . 



5.2 Gyroresonant surfing acceleration 

iKuramitsu and Krasnoselskikh! (|2005ah proposed a new type of surfing acceleration of 
ions resonating with a left-hand circularly polarized wave through the fundamental 
cyclotron resonance conditiorQ with n = +1, 

w - k x V R = n ci (14) 

It is known that a particle trapped by a monochromatic circularly polarized wave 
oscillates around the resonant point in a phase space. During a trapped oscillation, a 



4 In this and next subsections we consider the resonant interaction of ions and set the sign 
of u> positive for the left-hand polarized waves. 
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Fig. 8 From Iscnbcrg and Vasqucz (2007). (a) The cyclotron resonance condition between 
protons and parallel-propagating ion cyclotron waves. The solid lines plot the wave dispersion 
relation, while the dashed lines show the resonance condition for two particular proton parallel 
velocities with solid squares indicating the resonance points with the waves, (b) Resonant 
surfaces for protons in velocity space, (v z — oj/fc) 2 + v 2 ^ = const. The cyclotron interaction 
causes protons to diffuse along these surfaces, but not across them, (c) and (d) Similar to (a) 
and (b) but for a minor ion species with Qi < Q p . 



particle energy as well as a relative phase between the particle and the wave is bounded. 
However, if an external force is imposed to keep the relative phase constant, the particle 
can continue to gain energy from the wave electric field. 

Such a forced resonance condition is realized, for instance, in the foreshock region 
of a quasi-parallel shock where an electrostatic potential acts as the external force. 
Figure shows completely resonant (solid lines) and near resonant (dashed lines) ion 
trajectories in v x — x (a) and v z — Vy (b) phase spaces in a given electrostatic potential 
whose spatial profile is represented with a dotted line in (a). The two particles are ac- 
celerated monotonically in perpendicular direction while keeping the parallel velocities 
(almost) constant at v x ~ Vr(= 0.5). In their test particle simulation all the resonant 
particles are eventually accelerated regardless of their initial gyro phases due to the 
effect of phase synchronization which is shown in Figure [7] (c) in which time evolution 
of v y of the resonant particles having a variety of initial gyro phases (solid lines) and 
E y component of the wave electric field (dotted line) are plotted. 
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Fig. 9 Figures from Matsukiyo and Hada (2009). (a) Two four- velocity trajectories (black 
and gray curves) of relativistic electrons which are efficiently accelerated, (b) Relativistic fun- 
damental cyclotron resonance conditions (1 1 (3 6 for positive (+) and negative (-) parallel wave 
numbers. 



5.3 Cyclotron resonant heating of solar wind heavy ions 

It has been known that heavy ions in the fast solar wind have faster flow velocity and 
hotter temperature (especially in the direction per pendicular to Bp) than the p roton 
population, the main population of the solar wind. Ilsenberg and Vascmezl ([20071 ) have 
proposed a second-order Fermi acceleration mechanism for the preferential perpendic- 
ular heating of heavy ions by counterstreaming left-hand circularly polarized waves. 
In Figure [SJ the case for protons ((a) and (b)) is contrasted with that for heavy ions 
((c) and (d)): The second-order Fermi process works rather slowly for protons. This is 
because the signs of the resonance velocities for parallel and anti-parallel propagating 
waves are different, and the acceleration of protons by these waves does not occur si- 
multaneously. On the other hand, minor heavy ions whose cyclotron frequency is 
smaller than Q p have multiple resonant interaction points if the condition \v z \ < vq 
is satisfied (Figure [8] (c)jf|- Therefore, heavy ions can be accelerated by multi waves 
simultaneously as indicated by the thick arrow in Figure [Hid) . 



5 The authors took the z axis along Bo. The limiting velocity vo for the multiple interactions 
corresponds to the condition tangential to the wave dispersion relation (Figure [He)). 
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5.4 Relativistic particle acceleration in developing turbulence 



IMatsukivo and Hadal (|2009l ) found a very efficient acceleration of relativistic particles 



in developing Alfven waves in an electron-positron pair plasma, where the wave number 
and frequency spectra of the waves are developing through parametric decay instabil- 
ities. Figure [9] (a) shows two examples of trajectories of electrons which are efficiently 
accelerated along the vertical line (uii = 0). The gray line denotes an accelerated 
electron in the course of the decay instability self-consistently reproduced in a one- 
dimensional full particle-in-cell simulation, while the black line is obtained from a test 
particle simulation in modeled electromagnetic fields. This acceleration is interpreted 
in the following scenario: In the developing turbulence, there exist waves propagating 
both directions parallel and anti-parallel to the background magnetic field. The fun- 
damental cyclotron resonance condition (n = +1) for relativistic electrons is V x = Vr 
where Vr satisfies 

w - k x V R = ^ (15) 
7 

With the four velocity (uu,u±) = (v x "//c, v±^/c), (fT5|) is rewritten afl 

Figure [9] (b) shows the relation between uu and ttj_ given by (|16|l . where the frequency 
uj/Qce = 0.11 and the parallel and anti-parallel wave numbers k x c/f2 C e = ±0.65 are 
chosen. It is remarkable that at u\\ = the two resonance curves overlap so that 
an electron can interact with two counterstreaming waves at the same time and be 
efficiently accelerated. This scenario is similar to the one for heavy ions f £l5.3[) except 
that the former is made possible by the relativistic reduction of the cyclotron frequency 
(i? ce /7 < Oce) instead of the intrinsic difference in cyclotron frequencies for the latter 
(Qi < O p ). 



5.5 Relativistic turning acceleration 



Recently, it is suggested that relativistic electrons in the earth's radiation belt are 
coherently accelerated by large amplitude whistler waves satisfy i ng the relativ i stic cy - 
clotron resonance condition (|15[) (e.g., see lOmura et al.1 (|2007T l: iFuruva et all Hooi); 
iTao and Bortnik Large amp l itude whistler wave s are actually observed as 'cho- 

rus' emissions ( Cattell et al.l d2008T ); Cully et all d2008T)) . The essence of the 'turning 
acceleration' model ( Omura et al.l I 2007t l: Furuva et al.l ( 20081 )) is the relativistic re- 
duction of the cyclotron frequency in 1)150 . Figure flOl shows two cases with Vr < 
(with 7 = 7^;) and Vr > (with 7 = jd > 7s). An electron with 7^ and Vr < 
is assumed to start interacting with the (counter-streaming) chorus wave. Since the 
chorus wave has a large amplitude, the electron is trapped and accelerated to increase 
7 (i.e., reduce Vr) coming to the turning point with Vr = 0. After the turning point, 
the electron moves toward the same direction of the wave propagation, so that the 



Note the following relation for the Lorentz factor, 

7s{ l_(^ +U i)/ C 2 } -V2 =(1+ „2 +ui)1 /2, 
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Fig. 10 Figures from [Omura ct al. (200^). (a) Fundamental relativistic cyclotron condition 
1151 1. (b) Illustration of 'relativistic turning acceleration'. Large amplitude chorus waves are 
assumed to be excited near the equator and propagate from there to the high latitude region 
along the field line. An electron coming from the low altitude mirror point toward the equator 
starts interacting with the chorus wave. 



wave-particle interaction to continue long enough to further accelerate the electron up 
to 7 = 7£>, where the electron is detrapped. 



6 Concluding remarks 

We have given an overview on the importance of resonant wave-particle interaction 
processes in space plasmas. Although each of elementary processes are well known 
one by one, there is a necessity to unify these processes. For example, as we see in 
i|4] to consider interaction process between particles and oblique propagating waves 
we need to take into account of all resonance interactions with harmonic numbers 
n—Q, ±1, ±2, and sometimes even with subharmonic numbers n=±l/2, ±3/2, .... 
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In §5.1-3 we see the various roles of cyclotron resonant condition in the processes, 
shock injection, gyroresonant surfing, and ion heating in the solar wind. Of course, the 
considerations for the elementary resonance processes alone do not suffice for compre- 
hensive understandings of targeting physical processes: Injected electrons at the shock 
front ( i)5.1|l or surfing-accelerated ions at the shock front ( i)5.2|l should be accelerated 
further in diffusive shock acceleration process, for example, to acquire relativistic en- 
ergy. For the dynamics of heavy ions in the solar wind, the second-order Fermi heating 
mechanism ( >i5.3p is incorporated into an inhomogeneous coronal hole model with the 
effects of gravity, ch arge-separation electric field, and mirroring force in the decreasing 



solar magnetic field l lsenberg and Vasaued 1 200' 



In >i5.4l and M5.5I We further see that the relativistic resonance condition makes the 
interaction process very rich. Note that these considerations are for the behavior of rela- 
tivistic particles in nonrelativistic turbulences. While there are a few pioneering studies 
in stochastic acceleration proce ss in relativistic turbulences (e.g. JViratanen and Vainio.1 
( 20051 ) ; IO 'Sullivan et all |2009l )), further works on relativistic turbulence are needed to 



get comprehensive understan ding, f o r exa mple, of the acceleration process of ultra- 
high-energy cosmic rays (e.s.. Illillasl |l984|)). 

Since we have limited the discussion to the waves in cold plasmas, wave modes 
raising from the finite temperature effect are outside of our scope. Among these modes, 
kinetic Alfven waves have been a focus of extensive inv estigation. Here we r e fer th e 
interested readers to a few representative pa pers (e.g., Illasegawa and Chenl ( 197d ): 



Leamon et all l| 19991 ) : iLvsak and Lotkd l|l996r i 
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Appendix A 



For the j species (j=i or e), the plasma frequency is defined as uj p j = y/Annoe 2 /rrij , 
and the cyclotron frequency as Q c j — eBo/rrijC. The Alfven velocity Va is defined as 
Bo/\/47rno(mj +m e ), which is chosen to be 10 _4 c in the test particle calculation. The 
ion mass is taken 1836 times m e . 

The wave electric field E w oc_exp_i(kr — uii) should satisfy the following matrix 
relation (see, e.g., Chapter 2 of lstixl (|l992l )). 



ME W = 0, (17) 

with a matrix M given as 



P — N sin 6 Af 2 sin6»cos^ 
M = { S — N 2 -iD ) (18) 

7V 2 sin0cos0 iD S-N 2 cos 2 t 
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where N = kc/uu is a refractive index, and P, S, and D are defined as 

2 9 

P = P{ U ) ei44 

S=\{R{u) + L{u)) 

D=\{R{u)-L(u)) (19) 
R(u) = 1 - ^ e 



w(w + f2 c i) w ( w — ^ce) 
2 2 



Note that the condition detAf = gives the wave dispersion relation. 

From the first and second column of the matrix equation (|17|) . we have 

TV 2 sin 9 cos 8 
Ew ' X ~ ~P-7V2sin 2 e ^ 

tu ' 1/ - S - N 2 



(20) 



Now switching from the complex formulation to the real formulation, we write the xyz 
components of the electric field, (E w ,x, E w ,y, Ew,z), as 

(E w . x , E w . y , E w . z ) = {g x costp,g y sinip,g z costp) (21) 

where the relative amplitudes of (g x , g-y, g z ) are determined by (|20[) . and tp is the phase 
angle, 

ifi = k ■ r — ujt + ifo (22) 
where ipo is the initial phase angle, respectively. From the Faraday's law, we have 

B w = N (cos 9, 0, sin 8) x E lu (23) 

from which the xyz components of the magnetic field, (B w ,x> B w ,y, B WjZ ), are obtained 
as 

B WtX = — N sin 9E w .y = — Ng y sin9 sinp, 

B w ,y — N(sin9E WiX — cos9E w . z ) — N(g x sin 9 — g z cos 9) cos <p, (24) 
B WyZ = N cos 9E WiV — N g y cos 9 sin (p 
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Fig. 11 Figure from Kuramitsu and Krasnosclskikh (2005b) showing the definition of angles, 
rfip, <p w , and ip. 



Appendix B 

In a finite amplitude monochromatic electromagnetic wave propagating parallel to the 
background magnetic field Bp charged particles can be phas e -space tra p ped a r ound the 



resonant velocity (see e .g., Palmadesso and Schmidt! dl971 }j_ Helliwelll (1 19741): SchmittJ 



197d~):lKarpmanl ( 19741 ) jMatsumotol (|l979h ~ Hosmno and Terasawal Jl985h:lKuramitsu and Krasnoselskikhl 



2005b)). Here we follow the description in a recent article bv lKuramitsu and Krasnoselskikh] 
2005b ). Firstly, the gyrophase angle for an ion is defined as 4> p = tan^ 1 (v z /v y ), the 
wave phase angle as <j>w = kx — Lot + a w (a w : the initial value) , and their difference 
as ip = <f> w — 4>p (Figure [TTJ . Next, with the velocity of ions in the wave rest frame 
(u x ,u±) = {vx — oj/k,vj_) the cosine of the pitch angle is denned as /i = u x /\u\, where 
\u\ = (it^ + it^) 1 / 2 is the constant of motion (namely, the ion energy is conserved in 
the wave rest frame). The second integral of the ion motion \ can be written as, 

X= \ ( M + I) 2+6(1 -V 2 ) 1/2 ™si> (25) 

where k = ku/il c i, and b the wave amplitude normalized by |Bo|. With x, the equation 
of motion is simplified as 

Figure IT51 shows the ion trajectories in the [i — tp phase space with k = 3 and b — 0.01, 
0.1, 1, 10. When the wave amplitude is small ((a): 6 = 0.01), ions are trapped around 
(/i, ■;/>) = ( — 1/k, 7r), which is due to the linear cyclotron resonant interaction. As the 
wave amplitude b becomes larger ((c) and (d)), the trapping region also becomes larger 
and a new trapping region is born in a location far from the original point ( — 1/k, it). 
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Fig. 12 Figure from Kuramitsu and Krasnosclskikh (2005b) showing particle trajectories in 
the y, — ip phase space with n = 3 when (a) b = 0.01, (b) b = 0.1, (c) 6=1, and (d) b = 10. 
One curve corresponds to one constant of motion \. 
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